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1. Introduction 


Presented here are various flutter-analysis techniques ob- 
tained by approximating the aerodynamic-influence-coefficient matrix 
(within the flutter range of the reduced frequency) in a form which 
is convenient to write the dynamic equations of the aircraft in a 
State Space format i.e., of the type 

sx “ ^ (1) 

where s is the complex frequency, x is the state vector, A is a 
f" 'uency-independent matrix or may be treated as such (for instance 
if the problem is solved by iteration with A constant in each 
iteration) . 

In recent years optimal ttructural design and active control 
technology have gained an increasingly important role in the 
improvement of aircraft performance; sophisticated techniques are 
required, especially for flexible aircrafts. The optimal design 
as well as the control-systems design involve the interaction of 
numerous disciplines including structures, aerodynamics, aero- 
elasticity, guidance and control. 

It may be noted that for the analysis of flexible aircrafts 
the equations of the aircraft dynamics may be cast in the format 
of Eq. (1) (with constant A) by approximating the aerodynamic 

matrix with the two lowest terms of their Taylor series around 
k=0 (where k is the reduced frequency) : This approximation is 

used for instance in FCAP (Flight Control Analysis Program, Refs. 

1 and 2) which is a general purpose program for the analy.sis of 
flexible aircrafts of arbitrary shape with active control. The 
low frequency formulation however is too restrictive for analysis 
of flutter or even control-surface dynamics. In this case fully- 
unsteady-aerodynamic influence coefficients are being used in 
FCAP (Pefs. 1 and 2) . 


In addition it should be noted that there exist already 
aerodynamic codes capable of dealing with transient response 
(analysis off the imaginary axis) such as SOUSSA (Kefs. 3 to 8) 
which is the most general one presently available, since it 
can be used for steady, oscillatory and fully unsteady, sub- 
sonic and supersonic flows around aircrafts having arbitrary 
shapes. Note that SOUSSA is the only fully unsteady code for 
complex configurations.* 

However, the currently, available flutter- analysis routine 
are limited to the analysis of the aerodynamic loads for purely 
imaginary reduced f requeue '.es (oscillatory flow) . This short- 
coming is partly duo to the fact that aerodynamic loads ob- 
tained from existing aerodynamics calculation methods had been 
limited until very recently to sinusoidal motions only and 
hence the classical flutter analyses (such as the V»g method) 
are limited to the imaginary axis. It is obvious that the 
flutter analysis needs to be re-formulated if one wants to 
take full advantage of the transient response feature of the 
aerodynamic routine SOUSSA. 

The purpose of this paper is to present and assess 
various methods for approximating the aerodynamic forces 
so that the two basic features introduced above (State Space 
formulation and of f-the-imaginary axis analysis) are retained. 
The advantages of retaining these features are considerable, not 
only in simplifying the flutter analysis, but especially for 
more advanced applications such as optimal design of active 
control in \/hich the flutter is merely a contraint to the 
optimization problem. 


*SOUSSA is used as aerodynamic subfirograi' iaE’CAP.^'^ 


Therefore presented here are 3 procedures for the calcula- 
tion of flutter frequencies and speeds. In addition to the varia- 
tions of the traditional V-g method, a new formulation incorporat- 
ing state space techniques of the type x-^ is introduced. The 
advantages of the state space format are also presented. 
Preliminary results are briefly described. 
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2. V-q Method 

The V-g method is briefly outlined here. Consider the 

t ut C 

Lagrange equation of motion for simple harmonic notion, e ^ 

- lO” K *■ 0 t J \ * 

where x is the vector of the’ generalized Lagrangian coordin- 
ates, M and K are the mass and stiffness matrices, q is the 
dynamic pressure and E(k) is the aerodynamic influence matrix 
(function of the reduced frequency k“Wl/v)and g a fictitious 
damping coefficient introduced to avoid the use of complex 
frequencies. With some algebraic manipulations, a standard 
eigenvalue problem is resulted which has the form (Ref. 9) 

(A-ni) x=0 

where A is a function of k. it should be noted that n has 
the expression 

For different values of k, a corresponding fl results. 

The flutter frequency, hov/evor, is the one in which g=0, i.o. 
n is real. This yields w and hence V (from k) . 

Preliminary results arc presented in Figs, a and b and 
compared with the ones of Ref. 10. 
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3 . Truncated Power Seriog Approximation to the Aerodynamic 

IntTuonco Coefficients 

The flutter nr \ ’s can be greatly simplified (by cast- 
ing it into a space- ;e format) if the aerodynamic influence 
coefficients are expressed as a truncated power series in 
the reduced frequency 

E (l<) - ^ - I k ^ + CC < 2 .) 

Preliminary results for the approximation of the aerodynamic 
forces with three-terms power series are presented in Figs. 1 
to 18, which present the real and imaginary parts of the nine 
aerodynamic coefficients relative to the f^rst two bending 
modes (i«l and 2 respectively) and one torsion mode (i»3) for 
a rectangular uniform wing.* Note that in Fig. 1, and Fig. ^ 
the third order approximations replicate the actual valuei quite 
accurately. Unfortunately, not all the coefficients were 
represented precisely by the third order power scries: see in 

particular Fig. 5. In these particular cases, either higher 
order povar series or an iterative scheme must be used if 
accurate approximations are desired. 


* 




These values were obtained with the procedure described 
in Appendix A. 
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Natural State Space Method 


Combining Eqs. (1) and (2) one obtains, in the time domain 

M X ■» k X t 1 / — £ , r » 1 1 ^ i T ^ 

— V"--' V-* ' 

which can be recast in the State Space format as 



By varying the values of V, standard root-locus techniques 
can be employed to determine the roots as a function of the 
velocity, and therefore the stability boundary (flutter speed) . 

Unfortunately, the validity of the State Space method is 
contingent on the invertibility of the matrix. If is 
ill conditioned, then small truncation or roundoff errors will 
yield large variations in its inverse, modifying the dynamics 
represented by Eq. (2) . 


Proliminary results were not satisf Actory . The problem 
was traced down to the near-invertibility of the matrix. 
For a typical evaluation of the coefficient matrices of the 
power series, took on the following value 


8.15214 


-6.60933 

10.8949 

-7.32938 

• 

6.28561 

10.11678 

-31.13977 


24.86876 

39.64769 


Note that the value of the determinant is D«.08. Note also 
that the product of the three diagonal terms is approximately 
equ.?l to 2000. This implies a very strong elimination of 
significant figure, which is to bo expected, since the first row 
is approximately equal to the second multiplied by -1.04. 


5. Pade'B Space«»5tatc Method 

The rcprosentation of _K(s) can be cxproascd as a matrix 
Fade approximant in the form 

j E^‘>) J- ' [ft ♦ 1. B ♦ 1 

which could yield a better approximation at high frequencies 

than the power series approximation. 

# 

Introducing a new state variable: 

r ■- bi- ' [ 8- * * ’i' &. ] ; 


and modifying Eq. (5)^ results in the following state space 
model : 



Hntrices, A, B^^, are all dependent on If 

is near singular, the accuracy of the Fade approximation 
will be sensitive to truncation and roundoff error. The 
preliminary result. s were not satisfactory due, again, to 
the near singularilry of the the matrix. 


* 


4 

Neglecting terms of order k one obtains: 
A=e"^E 2 , Bq=-AEq, Bj^=Eq-AEj^, B 2=E]^-AE2 
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6. Iterative State Space Method 

In an effort to alleviate the problem elicited by this 
ill conditioned matrix, an iterative State Space method is 
considered. In order to justify the method, consider the 
simple equation 


f(t) {£ y 

The amplitude of one of the roots of f(x) approaches 
infinity as £ goes to zero. Hence, the root does not play 
an important role in determining the locus of roots obtained 
by varying the co'.'I ficients of the lower order terns^c(as 
long as it is in the left hand side of the plane) . 

The other two roots may be obtained by solving by 
iteration the quadratic equation 


or 


cr 


c\ fa ft fC = O 


■* !o < \ C =0 

- -2-3 

v;ith a»a + Ex, b*b + E x c * c+E x respectively. Therefore 
Eq. (3) is equivalent to (in the frev]uency domain) 

( M t E, ) X ♦ s fc. , X 4 (Cl 4 L^) A i <2 


X A 2 *** 

(v/ith + s E^) or, in space-state format 


H ' 1 - 1 

1 j 

1 o 

' f o -ill'. 



A 


( ». - 


O K ♦ 1 

b'S. £, 1(5.) 



v/here Xj^ = *ix. 


I 



It should be noted that the roots of Eq. (3) which are 
not obtainc ' are expected to be highly damped and therefore 

f 

of little interest. 

Also it should be noted that a variation of the above 
method is obtained if the actual value of E is used (by 
employing obvious appropriate definitions for Eq and in 
each Iteration) . 

Results using this method are now being obtained. 
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7. Concluding Romarka 

Several n^thods to solve the flutter equations have been 
briefly described. Preliminary results have been presented. 

The state* space method has certain advantages over 
V-g method. Because the approximation to the aerodynamic 
influence coefficients is a function of the complex variable 
s, the truncated power series can be evaluated at points off 
the imaginary axis. Hence, eigenvalues obtained which are 
not purely imaginary have physical interpretation. In the 
V-g method, eigenvalues obtained for which the value of g is 
not identically zero have no physical moaning, but are merely 
mathematically introduced. 

■•■n the V-g method, the value of k for a flutter crossing 
is .xiLSt determined, then V is computed. If it is desired to 
obtain the lowest value of V for which a crossing of the 
imaginary axis occurs, k must sweep out to infinity, for the 
redu'^ed frequency varies as 1/V. In the new method, V is 
directly inputed as a parameter, hence tiio sweep can be deter- 
mined at the first crossing. In the V-g method, many flutter 
crossings will occur. 

In addition, system theoretic concepts can be employed in 
the formulation of an active control solution to flutter. For 
example, optimal control techniques can be carried out subject 
to the dynamical constraints of the system. Also, expression 
of the aerodynamic influence coefficients in the state variable 
format enables the analysis of nonlinear dynamics using various 
methods such as multiple time scaling or the Lie transform. 
Dynamic stability of tl^e wing under study can be evaluated as 
v/oll as overall system performance. 


APPENDIX A 


EVALUATION OF IIATRICES E. 

_k 

Fig. 1 through Fig. 18 compare the polynomial approxi- 
mations with the actual values of the aerodynamic influence 
coefficients obtained by SOUSSA for each element of the force 
matrix. The value of the coefficients are plotted as a function 
of the reduced frequency, k. The squares represent the actual 
values^ while the triangles represent the polynomial approxi- 
mation . 

The coefficient matrices Eq and were obtained as 

E2 and were computed as 



The expression in brackets is constant if _E(k) was equal to 
a cubic polynomial within the range of frequencies (i»l,...,N) 


RHPRODUCIBILITY OF THE 
ORIdlNAL PAGE IS PCX>H 
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Flutter speed os a function of for a rectangular wing v/it'n 
AR = 16, M = 0, T = 0. 1%, n = 5, >'ju = 0.2, “ 0.5, and 

fvJX = 8 , NY =* 10. Results ere corrpared v/ith exact solution 
given by tv/o dimensional airfoil theory (Rof.io) (^ta “ “0.2C, 


two dimensional 


0 present method 


b. Flutter spoerl ns o function of ‘^\/’S>rj ^ rectangular v/Ing v/ilh 

AR ~ 16, h\ — 0, T = 0 . I*t;o, - 10 , X:v — 0.2, -0.5 or.d 

MX ■“ 8, NY - 10. Results ciu compa f<! v/Ilh exact solution 
given hv two dirrensionci airroil theory (Ref. 10)(Xr a = -0.2C). 
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